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Abstract. Static self-consistent methods usually allow to determine the most 
probable fission fragments mass asymmetry. We have applied random neck 
rupture mechanism to the nuclei in the configuration at the end of fission 
paths. Fission fragment mass distributions have been deduced from the pre¬ 
scission nuclear density distribution obtained from the self-consistent calculations. 
Potential energy surfaces as well as nuclear shapes have been calculated in the 
fully microscopic theory, namely the constrained Hartree-Fock-Bogolubov model 
with the effective Gogny DIS density-dependent interaction. The method has 
been applied for analysis of fission of 256,258pj^^ 252 q£ ^®®Hg and compared 

with the experimental data. 


PACS numbers: 21.65-l-f, 24.75.-|-i, 25.85.Ca 
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1. Introduction 

Fission is one of the dominant decay channels of the 
heaviest nuclei. Physics of this process is crucial 
in determining stability of heavy and super-heavy 
isotopes. One of the basic observables of fission, 
obtained directly in experiment, is fragment mass 
distribution. Measured yields allow to determine 
type of fission as well as to deduce the properties 
of the mother nucleus. Prediction of the charge, 
mass and total kinetic energy (TKE) distributions of 
hssion fragments is still a challenging task for the 
theory of nuclear fission. A number of attempts have 
been made to describe mechanism of sharing nucleons 
between hssion fragments since the spontaneous hssion 
phenomenon was discovered. 

Historically, in the hrst theoretical description 
of hssion the macroscopic liquid drop model was 
used. In this approach the competition between 
Coulomb repulsion and surface energy of deformed 
charged liquid drop of nuclear matter is analyzed. The 
mass distributions calculated within this method are 
symmetric, as a consequence of ignoring microscopic 
effects, that are responsible for octupole deformation of 
mother nucleus and deformations of the fragments [I]. 
The more sophisticated approach, scission point 
model [5], assumes that hssion properties can be 
derived from the scission point conhguration that is 
well dehned in the evolution of hssioning nucleus. The 
scission point is described as a conhguration of two 
deformed touching fragments. The potential energy 
is calculated using macroscopic-microscopic model. 
The probability of certain fragmentation is in inverse 
proportion to the scission point potential energy. This 
method allows to determine the probability, that the 
certain number of nucleons will be incorporated to 
one of the fragments from the energy of the scission 
conhguration. 

Within improved scission point model (HI Ql C] 
one can obtain distribution of hssion fragments and 
the mean value of total kinetic energy. In comparison 
to Ref. [2] this model allows to calculate energy of 
interaction between fragments and deformation energy 
at the scission point conhguration, what results on 
obtained mass and TKE distributions. The main 
disadvantage of this treatment is large number of 
phenomenological parameters, e.g. touching distance 
of hssioning system, surface tension coefficient with 
parameters htted to the magic nuclei. 

The authors of the microscopic scission point 
method [6] also deduce hssion properties as fragment 
mass and TKE from the analysis of dinuclear system 
that may be created after scission. Energies of 
nuclei are obtained in the microscopic self-consistent 
calculations. Depending on mass asymmetry and 
deformation of the fragments energy of the system is 


calculated and hssion probability is determined. In 
this approach strong assumptions on distance between 
fragments and deformation of the fragments at scission 
point have to be made. 

Another approach includes dynamics of hssion in 
the analysis. It treats the nuclear shape evolution 
as a Brownian motions of nuclear system on the 
potential energy surface (PES) [HIH]. The direction of 
the motion in the hve-dimensional deformation space 
is determined randomly using Metropolis’ method. 
The considered PES is calculated within macroscopic- 
microscopic model, where the potential energy consists 
of liquid drop part (with deformation-dependent 
coefficients) and microscopic corrections. During 
“random walk” on this surface each shape might be 
obtained with accordance to its statistical weight. 
The dimensional depletion causes that agreement 
with experimental mass distribution is getting worse. 
The additional parameter - critical neck radius was 
introduced, and its value results much on the obtained 
fragment mass yields. 

There were also several attempts to describe 
fragment mass distribution in a fully microscopic 
way, i.e. using Adiabatic Time-Dependent Hartree- 
Fock-Bogolubov (ATDHFB) method [5] or Hartree- 
Fock-Bogolubov (HFB) method with Skyrme Energy 
Functional [10]. ATDHFB model allows to include 
dynamical effects to the description of the fission 
process. The probability of mass division is 
proportional to the total flux of the wave function 
through the scission point in a particular configuration. 
Including dynamic effects gives much better agreement 
with data in comparison with static calculations. 
Especially the broadness of fragment mass distribution 
is closer to the experimental one. 

The HFB method with Skyrme energy functional 
was applied for description of induced fission process 
[10] . The authors discussed the impact of the triaxial- 
ity and discontinuity of the PES in the scission point re¬ 
gion on the fragment mass distribution. Quasi-particle 
occupations in the nascent fragments were localized 
and attributed to each fragment. However, no com¬ 
parison of theoretical and experimental distributions 
were performed. 

Dynamic calculations at low excitation energy 
were performed using Langevin approach m- Shell 
and pairing effects, as well as dissipation and 
fluctuation were included in the model. Out of 
fragment mass distribution it was also possible to 
obtain the time scale of fission process. The 
authors made a strong assumption, that both nascent 
fragments have the same deformation, what argues 
with several observations of shapes of fission products 

m- 

Published lastly GEF code [13] predicts surpris- 
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ingly good fission fragment mass yields with a very 
simple construction of the PES of fissioning nucleus. 
Namely, macroscopic potential is modified by parabolic 
corrections to simulate shell effects. Once fitted it re¬ 
produces wide range of experimental fission results. 

The main inspiration for our investigations is the 
idea proposed by U. Brosa et al [H HSl Hi- The 
authors assume direct dependence of fragment mass 
distribution on the random neck rupture mechanism of 
mother nucleus in its pre-scission shape. The method 
describes splitting of a nucleus as a consequence 
of hydrodynamic fluctuations induced by random 
vibrations of nuclear surface. The pre-scission 
deformation of the fissioning nucleus has a decisive 
influence on the fragment mass yield. The probability, 
that the rupture of the neck occurs in a certain 
position, decreases with the radius of a neck. In 
this method no assumption on properties of fission 
fragments is needed. This simple macroscopic model, 
allows to describe many fission properties as the 
TKE distribution and the dependence of neutron 
multiplicity on fragment masses. 

We apply the HFB model, employing finite- 
range Gogny forces to calculate the PES and mass 
distribution during fission. The previous investigations 
performed with this model proven, that it describes 
with reasonable accuracy the experimental data of 
the most relevant properties of unstable nuclei for 
different decay modes, such as spontaneous fission 
[HI m m HO] and exotic process of cluster 
radioactivity [52] . Also the asymmetric fission of i®°Hg 
isotope was successfully explained within microscopic 
approach [23]. In this paper detailed analysis of 
the scission point configuration determined in the 
microscopic calculations is performed. From the single 
mass distribution we deduce possible fragment mass 
asymmetries. 

This paper is built as follows: In the second 
section details of calculations are described. Results 
obtained for four isotopes are presented in the third 
section. Finaly fourth section include our conclusions 
and discussion of the results. 

2. Method 

The analysis of the scission point configurations 
requires realistic description of the nuclear matter 
mass distributions at large deformations. They 
may be provided by the calculations performed 
within the self-consistent constrained HFB model 
with the effective Gogny density-dependent interaction 
[55] . The popular DIS parametrization [55] was 
used. The computer code of m was applied for 
numerical calculations. The PES was determined in 
the deformation space of quadrupole and octupole 


deformations. In the PES the fission paths were 
determined as continuous lines connecting local 
minima of energy for fixed value of quadrupole 
moment. In each point of the constrained calculations 
nuclear matter spatial distribution can be easily 
derrived from the single-nucleon wave functions. 

2.1. The pre-scission point configuration 

The scission point is defined as the configuration of 
nuclear system in which molecular shape (i.e. two 
fragments connected by a thin neck) of a mother 
nucleus converts into two separate fragments. In some 
simple models it is represented by the configuration 
of two touching nascent nuclei. In fact, when a 
realistic leptodermous density distribution is taken 
into account, one can find, that during the splitting 
nuclear matter density in the neck decreases gradually 
from the bulk value to zero while surfaces of the 
fragments overlap. Thus we can not define one scission 
configuration but rather a session region in which two 
fragments are being separated. 

In the self-consistent calculations of the PES it 
is easy to locate the ranges of the scission region. 
Increasing constraint on deformation parameter, e.g. 
quadrupole moment, leads to more elongated shape of 
a nucleus with a thinner neck. At some point sudden 
change take place. The energy minimization procedure 
leads not to molecular shape but a solution with two 
nuclei separated by a few fm distance is found. Such 
dinuclear system usually has got energy much lower 
than the compound nucleus calculated in the previous 
step. Surely, scission take place between these two 
configurations. In this way the last point on the fission 
valley before rupture should be called pre-scission point 
and the next point with two separated fragments post¬ 
scission point. The line connecting pre-scission points 
on the PES should be called scission line (or more 
precisely pre-sisedon line), see Fig. [T] Intermediate 
configurations of the scission region are hard to obtain 
in the self-consistent procedure as they require multiple 
constraints, including e.g. neck parameter [221123] . 

The evolution of shape along fission path goes 
smoothly with increase of deformation up to the the 
pre-scissin point. It was found [m El] that many 
fission properties are determined by the shape of the 
PES from saddle to scission, where a nucleus takes a 
molecular shape. Especially the configuration in the 
pre-scission point is crucial in the fission process as no 
further evolution of a compound nucleus is allowed. 

In most cases determination of the scission line 
is trivial, as in the constrained calculations many 
observables essentially change their values: e.g. energy, 
hexadecapole moment, neck thickness. Such line can 
be easily noticed in Fig. [IJbetween Q 2 = 130 b and 250 
b as rapid change of energy of the system take place 
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Figure 1. The PES of The white lines correspond to 

the fission paths. 

along scission line. 

In the case of symmetric fission path in ^^®Fm 
at Q 3 = it is not so clear as molecular shape 

converts quite smoothly into dinuclear system without 
sudden drop of energy. It can be easily found in 
the panel (a) of Fig. [2l where fission path smoothly 
converts into hyperbolic decrease characteristic for 
Coulomb repulsion between separated fragments. Thus 
fall of energy can not be a criterion to determine the 
scission point. Detailed calculations show only a tiny 
kink with change of energy derivative when the neck is 
about to disappear. To prove, that this is real scission 
point, first we calculate neck parameter Qjy defined as: 

Qn = J Piz,r±) exp (1) 

where ao is a parameter describing width of a neck 
region with an arbitrary chosen value oq = 1 fm. 
Qn describes number of particles in the neck region. 
Middle panel of Fig. [2] shows changes of a neck 
parameter along the fission path. In the first part of 
fission path Qn decreases slowly, whereas beyond Q 2 = 
129 b it drops down rapidly. It means that sudden 
change in density distribution take place here. Neck 
parameter takes small positive and almost constant 
values from <32 = 130 b. This is characteristic 
behaviour of Qn for two separated nuclei when 
only tails of density distribution of both fragments 
contribute to Qn- 

As it was mentioned, scission point may be defined 
as a configuration, where rapid change of matter 
distribution from molecular to dinuclear system occurs. 
To check discontinuity in changes of density of nucleons 
it is useful to follow procedure of Ref. We may 
define the density distance function Dpp, given by: 

Dpp. = J \pir) - p'{r)\dr, (2) 



120 122 124 126 128 130 132 
Q2 [b] 

Figure 2. The fragment of the compact fission path of ^^^Fm 
in the vicinity of scission point. (a) The energy, (b) the 
neck parameter Qjv s,nd (c) the density distance function Dpp/ 
between two neighbouring configurations as a function of the 
quadrupole moment Q 2 - 

where p(r) and p'(r) are the the local spatial densities, 
calculated in two following steps during the nuclear 
shape evolution. This function takes small values 
during smooth evolution of nuclear shape. When two 
different configurations are compared it is manifested 
by a sharp single peak. 

In the panel (c) of Fig. [3] we can find again that 
function Dpp, is four times larger at Q 2 = 130 b in 
comparison to the neighbouring deformations. We can 
deduce rapid change of configuration between solutions 
at Q 2 = 129 b and at Q 2 = 130 b, namely rupture of 
the neck is observed. 

All above mentioned arguments show that at Q 2 = 
129 b pre-scission point can be found on the symmetric 
fission path of ^^®Fm. 

2.2. Random neck rupture mechanism 

It was proven in Ref. [331 [33] that the structure and ba¬ 
sic properties of the nascent fragments are preliminary 
determined in the pre-scission configuration of fission- 
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ing nucleus. In the pre-scission point two fragments are 
already created in the form of two nuclei connected by 
a neck. The most of nucleons are localized in the frag¬ 
ments, while some of them (up to 20 ) still create the 
neck. They shall be incorporated to the nascent frag¬ 
ments during scission. Usually it is assumed, that the 
neck is ruptured in its thinnest place. In this way one 
can calculate how nucleons are shared and the most 
probable mass asymmetry can be found. Nevertheless, 
following the idea presented in Refs. [T51[T6], the prob¬ 
ability P of rupture of a neck depends on the energy 
Ecut needed to create a cut in the considered position 
z along the symmetry axis. For each position z we can 
calculate particular fragment mass asymmetry. 

The probability P of the rupture of a neck, leading 
to fragment mass asymmetry A 1 /A 2 reads: 

P{Ai,A 2 ) exp(-£'™t/r), (3) 

where T = \fV2M^EjA is temperature of the pre¬ 
scission deformation. Excitation energy at scission: 

E^^ = Eg.,.-EZf, (4) 

which was gained during the evolution from the ground ^ 
state to the pre-scission deformation. E^ut depends ^ 
on the linear density (t(z) of a neck in a position 
of z that corresponds to required fragment mass 
asymmetry. Note that two values of z give the same 
mass asymmetry but with lighter and heavier fragment 
placed at opposite sides. Cut energy is defined as: 

Ecut{z) = 2ja{z), (5) 

The expression for surface tension coefficient 7 = 
0.9517[1 —1.7826(1 —2Z/A)^] was taken in a from given 
in Ref. [53]. The cross section of a neck is given by: 

poo 

cr(z) = 27r / r_Lp(z, r_L)dr_L , (6) 

Jo 

where p(z,rj_) is the spatial density of nuclear 
matter. Finally, the probability for the rupture, with 
corresponding division of nucleons between fragments, 
takes the form: 

A 2 ) = exp[- 27 CT(z)/T] (7) 

This method is visualized in Fig. [3] where the 
case of ^®®Fm symmetric mode is considered. The 
top panel shows the shape of the nucleus at the 
half-bulk density p = 0.08 fm“^ in the pre-scission 
point HZ]. Each neck’s rupture position at the z 
axis is correlated with specific mass division between 
fragments shown in the middle panel of Fig. |3J The 
probability of the rupture decreases strongly, while 
the neck becomes thicker, according to Eq. (jT]). In 
consequence mass yield presented in the lower panels of 
Fig. [3] is quite narrow with the most probable position 
of a rupture corresponding to the equal split into two 
^^®Sn isotopes. 


Qp=129b 

Q3=0 



-10 -5 0 5 10 

z[fm] 



A1 


Figure 3. The upper panel shows the shape of the 
at a density of po = 0.08 fm~^ in the pre-scission point on 
the symmetric path. The number of nucleons as a function of 
symmetry axis 2 (middle panel) and the probability of division 
of nucleons between nascent fragments, depending on the cross 
section of the neck in a certain position calculated using Eq. 0 
(lower panel). 

3. Results 

It was shown, that fission in ^®®Fm represents specific, 
bimodal character [301 EH. Its mass yield is very 
narrow with a single peak. The TKE distribution 
is compounded of the low and high energy modes 
with equal abundance. They can be linked to the 
theoretically described asymmetric and symmetric 
fission channels on the PES, respectively m, which 
are visible on the PES in Fig. |TJ A nucleus in its 
ground state is quadrupole deformed {Q20 = 16 b). 
The symmetric path goes along reflection symmetric 
shapes and terminates at deformation close to Q20 = 
129 b, as it was discussed in the previous Section. The 
density distribution corresponding to this pre-scission 
configuration is shown in the upper part of Fig. [4^. 
The shape of a nucleus can be described as a molecular 
system of two spherical ^^"‘Cd isotopes connected by a 
thin neck containing 4 protons and 6 neutrons. The 
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fragments conserve N/Z ratio of mother nucleus and 
reproduce the shape and mass of the outer part of 
the fragments before scission. In panels (c) and (d) 
of Fig. U density profiles of ^®®Fm and two ^^"^Cd are 
compared along z axis and along perpendicular axis at 
the thickest part of the fragments. Out of the neck 
region very good agreement can be noticed. 





z[fm] 



a [fm] 


Figure 4. (a) Density distribution of in configuration 

close to the scission on the symmetric fission path (Q 20 = 129 
b, Q 30 = 0 b^/^) compared with the density distribution of 
^^■^Cd (spherical shape), (b) Linear density of ^^®Fm in its 
pre-scission shape (continuous line) and ^^“^Cd (dotted lines). 
Profiles of density for cross sections taken in (c) rj^ = 0 fm and 
(d) 2 : = —7.1 fm, 2 : = 7.1 fm. 





z[fm] 



r^[fm] 

Figure 5. Similar as in Fig. [4|but for pre-scission configuration 
on asymmetric path (Q 20 = 238 b, Q 30 = 66.7 b^^^). 
Density distribution and profiles are compared with ^^^Sn (in 
its spherical ground state) and ^®^Mo (Q 20 = 6 b). 

The alternative asymmetric path arises behind the 
barrier from Q 2 Q Ri 90 b and Qao sa 20 b^/^ and leads 
towards configuration close to Q 20 = 239 b, Q 30 = 
67 b^/^. The analysis of density profile of asymmetric 
pre-scission deformation is shown in Fig. [SK- This 
molecular configuration consists of spherical double¬ 
magic ^^^Sn and prolate ^°®Mo. In this case, a neck 
contains 8 protons and 12 neutrons. Again, in panels 
(c) and (d) of Fig. [5l nice reproduction of the density 
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Figure 6. Fragment mass distribution for the spontaneous 
fission of ^^®Fm obtained within presented method. The sum 
of distributions (dotted line), corresponding to the symmetric 
and asymmetric paths is compared with the experimental yield, 
taken from Ref. m- 

profiles of mother nucleus by nascent fragments can be 
found. 

In the symmetric mode 10 neck nucleons should be 
shared between fragments in the scission point. The 
most probable split of symmetric pre-scission shape 
leads to the production of two ^^®Sn isotopes. In the 
asymmetric fission channel ^"^^Cs and ^^®Rh isotopes 
are most likely created. If we assume, that neck 
does not have to be cut in the thinnest place, these 
numbers would convert into stretched fragment mass 
distribution. To this end the procedure presented 
in the previous section is applied twice, both for 
symmetric and asymmetric pre-scission shapes. From 
the density distribution linear density is calculated in 
each case, what can be seen in Figs. IHd and[5]D. Next, 
fission probability as a function of mass asymmetry 
is calculated from Eq. ©• The results are combined 
with assumption of equal abundance of both modes 
and plotted in Fig. |6] by dotted line. In the same Figure 
experimental results taken from Ref. |31j are shown. 
One can find, that symmetric peak is much more 
narrow than the experimental one. Asymmetric mode 
produce two additional side peaks at Ah ~ 143, k, 
115 which are not visible in the experimental data. 
Nevertheless asymmetric peaks fit in the region of the 
tail of mass yield of ^^®Fm. 

The main reason of inconsistency is lack of 
dynamic effects in the present analysis. In the final 
part of evolution of a nucleus, the PES is quite soft 
in Qa direction and wide fission valleys are developed 
(see Eig. [T|). We have considered only two pre-scission 
configuration at the end of fission paths. However, 
there is a big chance, that nucleus evolves to the 
other pre-scission configurations. Exact probability 
of reaching each point along scission line on the PES 


should be obtained in dynamical calculations. From 
each pre-scission shape fragment mass distribution can 
be deduced. Combination of these two effects would 
lead to broader mass distribution. 




Figure 7. Similar as in Fig. |4]but for pre-scission configuration 
of 256pj,o (Q 20 = 230 b, Q 30 = 68b^A). Density distribution 
and profiles are compared with ^®'*Te (Q 20 = 1 b) and ^®®Mo 
{Q 20 = 6 b). 
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Figure 8. Fragment mass distribution for the spontaneous 
fission of ^®®Fm isotope in comparison to the experimental data, 
taken from Ref. [32]. 

The fission properties of ^^®Fm differ substantially 
from its neighbour isotope with two more neutrons. 
Its spontaneous fission half-life is seven orders of 
magnitude longer and symmetric fission channel is not 
active. The detailed explanation of these facts can 
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be found in Ref. m- In only the asymmetric 

channel of fission is visible in the experiment. Analysis 
of the density profile of the pre-scission configuration 
(Fig. ED indicates, that the nascent fragments are 
close to the ^^'*Te and ^°®Mo with 6 protons and 
8 neutrons in a neck region. From the linear 
density function (t(z) presented in Fig. Et* the 
fission fragment mass distribution are obtained within 
presented method. The results are shown in Fig. [5] 
together with experimental mass yields. The most 
probable mass of light fragment is reproduced with 
good accuracy. The heavy fragment peak is slightly 
shifted in comparison to the experimental one. The 
broadness of calculated distribution is reduced. 



-15 -10 -5 0 5 10 15 20 


z[fm] 


Figure 9. Similar as in Fig. |4]but for pre-scission configuration 
of ^®^Cf (Q 20 = 217b, Q 30 = 66.5b®/^). Density distribution 
and profiles are compared with ^®^Zr (Q 20 = 5.5 b) and ^®®Sn 
(spherical). 


In Fig. in the density profile of ^®^Cf in its pre¬ 
scission deformation is shown as well as its linear 
density function. The fragment mass distribution is 
presented in Fig.lTOl In this case, the most abundantly 
produced fission light fragment is not reproduced 
within presented method. The predicted mass 
distribution of the lighter fragment is shifted towards 
heavier masses in comparison to the experimental 
data. Probable explanation of these discrepancies 
comes from the fact, that post neutron emission 
measurements are presented [32] and average neutron 
multiplicity is larger in the region of the most probable 
lighter fragment [34]. Similar to the results of ^®®Fm, 
theoretical mass distribution is too narrow, but the 



80 100 120 140 160 


A 

Figure 10. The same as in Fig. m but for 252 isotope. 
Experimental data were taken from Ref. m- 


peak-to-valley ratio is in good agreement to the 
observed yield. 




Figure 11. Similar as in Fig. |4]but for pre-scission configuration 
of iSORg (Q 20 = 130 b, Qao = 30 b^A). Density distribution and 
profiles are compared with ®'^Zr (spherical) and ^^Ge (Q 20 = 
8 b). 


The recent experimental studies brought unex¬ 
pected observations of /3-delayed fission of ^®°Hg [3^ 
135] . The most abundantly produced fission fragments 
have got fragment mass asymmetry Ah/Al — 100/80 
instead of expected fragmentation leading to two magic 
®°Zr isotopes. Detailed microscopic analysis [23] ex¬ 
plained the dominant character of asymmetric fission 
valley. The ^®°Hg, as a product of electron capture 
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Figure 12. The same as in Fig. [8] but for isotope. 

Experimental data were taken from Ref. m- 

in is created with excitation energy not larger 

than 10.44 MeV. In this nucleus the pre-scission point 
is in the vicinity of a saddle around 12 MeV above the 
ground state. Scission point excitation energy should 
be taken equal maximal energy available for the nu¬ 
cleus in its ground state. Density distribution in the 
pre-scission point as well as function a{z) are plotted 
in Fig. [TTJ Fragment mass distribution is compared 
with experimental data in Fig. 1121 The most proba¬ 
ble heavy and light masses {Ah = 100 and Al = 80) 
are well reproduced. Calculated yield is too narrow in 
comparison with the observed one. The peak-to-valley 
ratio is not reproduced as well. 

4. Summary 

The presented results confirm, that microscopic de¬ 
scription of the pre-scission configuration provides a 
lot of information about physics of fission. Applying 
macroscopic method proposed by Brosa to the self- 
consistently calculated nuclear matter density distribu¬ 
tion, we have deduced ambiguous fragment mass yield. 
The major characteristics of experimentally measured 
mass distributions, i.e. the most probable masses of 
heavy and light fragment are reproduced with reason¬ 
able accuracy. 

The broadness of yields is too narrow as dynamics 
of the process was omitted in current investigation. 
The value of temperature used in our studies may 
also affect the results. Between pre- and post-scission 
configuration energy usually significantly decreases. In 
consequence one may assume larger values of excitation 
energy and temperature, which leads to broadening 
peaks of mass distribution. In the neck rupture 
mechanism we have also ignored quantal properties 
of nucleons, which may impact on fragment mass 
distribution. We have shown, that even in asymmetric 


fission one of the fragments is close to spherical double¬ 
magic nucleus. Strong shell effects in the fragments 
as well as quantal nature of the neck nucleons may 
modify macroscopic model of random neck rupture. 
All presented distributions, deduced from our static 
calculations, are much narrower in comparison with 
experimental ones. Further investigations require 
inclusion of dynamical effects. It would allow to 
reproduce required broadness of mass yields. 
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